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Neutrino Oscillations in Gravitational Fields 
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' We study the quantum mechanical phases associated with the propagation of neutrinos in 

, gravitational fields under the WKB approximation for Dirac equation. Applying the result, 

OA . the investigations are made for gravitational effects on the neutrino oscillations, the lensing 

» I ' of neutrinos and so on. 
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§1. Introduction 



The neutrino oscillation seems, nowadays, to be established experimentally as a 
fact in nature. The physics associated with the neutrino oscillation is still interesting 
^ ' subject to study from both theoretical and experimental points of view. Among 

. those, the gravitational effects on the oscillation are, sometimes, considered to lead 

insignificant corrections, since those are very weak for the present status of neutrino 
observation. Nevertheless, such an effect inspires several physical interests to us. 
CD . For example, the breakdown of the equivalence principle was often discussed 

from the viewpoints of the influence on the oscillation as a matter like effect and of a 
basic problem in gravity. In those approaches, the breakdown is usually introduced 
by assuming tiny mass dependence of gravitational coupling constant, though the 
, way of the breakdown is not unique and may not always be natural, since the break- 

^j' down is usually set after the standard geometrical approach to gravity that holds 

^ . under the equivalence principle. 

Another interesting gravitational effect on neutrinos may be the gravitational 
lensing or the gravitational redshift, which can be observed as real effects for pho- 
tonsj-^' in this case, the equivalence principle must be assumed to derive those effects. 
5^ I One can expect to observe a similar gravitational effect in neutrinos too.^^ Namely, 

neutrinos will serve as an observational tool in astrophysics; and conversely, astro- 
physics will give an insight into neutrino physics. For example, if it is possible to 
observe simultaneously the deflection of the neutrinos and that of photons started 
from the same astrophysical object, then the masses of neutrinos will be determined 
instead of those differences. 

To study the neutrino oscillations in gravitational fields,'^''^ we have to somehow 
evaluate the Dirac wave functions in curved spacetime. As a useful way to do this, 
the forirP ^ J"^ dx^Pf^{x) is sometimes used as the heuristic gravitationally induced 
wave function phase without saying on its validity. Here, p^{x) is the four-momentum 
of neutrino at x. The wave function with such a heuristic phase is obtained in an 
approximate solution of the Dirac equation in truth. The purpose of this paper 
is, thus, to formulate the propagation of neutrino in gravitational fields by using an 
approximation like the WKB method for the Dirac equation, which was discussed by 
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Pauli— to study the electrons in external electromagnetic fields, with the application 
to solvable examples. The heuristic gravity induced phase becomes to have clear 
meaning in this approach. In the present paper, the equivalence principle is also 
assumed to hold. 

The formulation will be done in the next section within the framework of two 
generations for sake of simplicity. In section 3, we apply the resultant wave functions 
to neutrinos in gravitational fields associated with the Schwarzschild exterior and 
interior solutions. Section 4 is devoted to discussions: the gravitational correction for 
oscillation phases, deflection angles of neutrinos and its relation to that of photons, 
and other discussion on practical problems. 

In appendix A, we give an introduction to the approximation of the Dirac equa- 
tion in a curved spacetime according to Pauli's method. Some calculations associated 
with the interior Schwarzschild geometry are also given in appendix B. 

§2. The phases of neutrinos propagating curved spacetime 

The neutrinos propagate vacuum oscillating among their flavor eigenstates l^a), 
(a = e,fi,T), when those are superposition of mass eigenstates {i = 1,2,3). 

Writing the states in momentum representation by |z>), the flavor eigenstates can 
be expressed as the superposition \Da) = Yli^ai\i^i) using a global unitary matrix 
[Uai]- Here, \Di) is the state satisfying mi)\Di) = and = Sij for pi = 

{E,p^), {pfscillat = E"^ — mf) *\ Then, the evolution of those states \va) from xa 
to in vacuum give rise to the phases so that 



\va{B)) = ^Uaie-'P-^''^--^^i),) (2-1) 



In what follows, for simplicity, we confine the discussion to the framework of two 
generations; and so, the unitary mixing matrix has the form 



[Ua 



COS 9 sin 6 
— sin 9 cos 9 



(2-2) 



The phases in Ea. H2-ip acquire another forms when neutrinos pass through matter 
field or a curved spacetime with a metric g^i, ^ ri^^- In the latter case, the Dirac 
equations for the mass eigenstates of neutrinos become 

[i7"ef^)D^-m,c]|z.,) = 0, (2-3) 

where e^^^ and = d^ — ^<^°^^'^{a)(0)ij. &re inverse vierbein and the covariant deriva- 
tive for local Lorentz transformations, respectively. 

Hereafter, we will discuss the case that the spin connections iO(a){f3)fi a-re negli- 
gible compared with the derivative term 3^ **\ Then, the |z/j) in Eq. (|2-3|) can be 

*' We use the unit li — 1 except Appendix A. 
**^ The ratio of to 9 is the order of ^ where Vg is the gravitational scale, such as 

the Schwarzschild radius, of an astrophysical object, to which neutrinos pass around with distance 
7?. Here, E is the energy of neutrinos; and for example, the ratio becomes <^ 10~^^ for E ~ WGev 
with R ~ 10^ km, the sun like radius. 
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solved at xb in the form \iyi{B)) ~ e'^'-^^^^lufiB)) in the sense of WKB approxi- 
mation (Appendix A). Here, Si{B,A) is a solution of the Hamilton- Jacobi equation 

g''''{d^S,){d,Si)-{m,cf = (2-4) 
at Xb', and \uf{B)) is a state satisfying 

{rel^d,S, + m,c)\iyf{B))=0 (2-5) 

at the same spacetime point xb- We, here, require Si{B, A) 0, {xb xa) so that 
Si{B,A) tends to —pi ■ {xf, — x^) in flat spacetime limit. In other words, {vii^B)) will 
tend to e~*^''^^^~^'^^|Pj) in that limit. Thus, in the curved spacetime, we can write 
the counterpart of Eq. H2-l|) as 

Wa{B,A)) c^Y.^a^e''^^'''^^W^{B)). (2-6) 

i 

The equality ~ means that |fa(i?,^)) is an approximate solution of the Dirac 
equation at xb disregarding hd^luf^B)) as a negligible quantity compared with 
dfj_Si{B, A). It should be noticed, however, that, if xb is a point in the asymp- 
totic region such as e^^^ = 6a and df^Si = const., we may identify \i^j'{B)) with 
\Di{B)), although 5j's are functions of xa and xb- In such a case, Ea. ()2-6() is a right 
solution of the Dirac equation in the curved spacetime. 

Hereafter, we assume that the xb is an observational point of neutrinos located 
at flat spacetime region and apply to {|f°(-B))} the same normalization as {\i'i{B))}. 
Therefore, we obtain the following expression for the transition probability in the 
curved spacetime: 

P(a^b) = \{hm'^aiB,A))\^ 

= (2-7) 

In particular, we have 

P(z.e ^ z^e) = 1 - sm\2e) sin^ (^^^ , {AS = Si - S2) (2-8) 

We, further, confine the discussion to the neutrinos go through a static gravitational 
field on the way. Then the energy of neutrinos conserve; that is, the time variable 
in the phase function Si can be separated at any spacetime point x = (t, x) as 

Si{x) = -Et + S,{x), (2-9) 

In this case, the AS = Si — S2 in Ea. ()2-8() is replaced by AS = 81 — 82, and the 
transition amplitude becomes time independent. 

§3. Gravitational effects in neutrino oscillations 

We, here, discuss the phase Si of the mass rrii neutrino in Ea. (|2-6p due to the 
static spherically symmetric gravitational field produced by a spherically symmetric 
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mass M body of radius R. We assume that the body is at rest against both the 
source of neutrinos at A and observer at B. Then we may take the body as the 
origin of spatial coordinates, the spherical coordinates {r,6,(j)), in which the A and 
B are located at places with the radial coordinates Ra and Rb, respectively. 

To solve the Si, we have to study separately the equations for Si, which hold 
outside and inside the body producing gravitational field. 



i) Si for the exterior Schwarzschild solution. 



The Schwarzschild metric, which holds outside the spherically symmetric body 
at rest, is defined by the diagonal vierbeirPI^ 



rsini?, e 



(3-1) 



where = r,x'^ = ^,x^ = (j) are spherical coordinates, and Vg = 2GM/c^ is the 
Schwarzschild radius of the body. Since, the angular momentum L of neutrinos 
conserve in the spacetime under consideration, we may choose z axis as the direction 
of the angular momentum. In other words, the neutrinos are on the plane with 
7? = ^. Then, we can rewrite the phase function S in Eq. (|2-9j) as 

Si = -Et + Li(l) + Sir{r). (3-2) 

We note that each neutrino Ui may have different angular momentum Lj, since the 
momenta of respective particles pi = \/{E/cY — i'm-ic)'^, (f = 1, 2) at xa are different 
for TTT-i 7^ 771-2. Substituting the expression H3-2p for Ea. (|2-4() . we have 

2 



1 



-1 



1 



dSir 

dr 



14 



ImiC) 



0. 



(3-3) 



This equation can be integrated to give 




^2 




(3-4) 



under the approximation up to the order of Here, + and — represent the integral 
along a r increasing path and a r decreasing path, respectively. Prom this expression 
for Sir, tfie lower bound of r integration becomes pi = Li/pi, which should be 
read as the impact parameter of neutrinos determined at xa] and so, we can put 
Pi = P2 = Po- Then writing r = r/po, fg = rg/po and so forth, the right-hand side of 
Ea. (|3-4() becomes 
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Fig. 1. The path from A (source) to B (observer). 



under the condition Ra,b ^ Po ^ ^g. Here the Ra and Rb are radial coordinates 
of neutrinos at initial and final points, respectively. The integral of each term in 
Ea. (|3-5() can be simply carried out, and we have 



k=A,B L 



1 — cos 



-1 



Rk. 



PiC 



+ cosh-(i?.)-?fji-4 



(3-6) 



Now, 



— 1 are respec- 



4°^ = = A,B) and Ub = ^Ek=A,B 

tively (j) coordinates of A, B and the distance between A, B in the flat spacetime. In 
a curved spacetime, the (j) coordinate of neutrino Vi as the function of r is defined 
||l = + ||r = 0. Then (t)i^AB = -^^^^tT^ becomes the difference 0b - (pA 
for the mass rrij neutrino, to which a little calculation leads to 



'H,AB 



I 



^AB 



dr 



+ 



Li 



Li. 
2r 



Pi 



2Li 



k=A,B 



+ 



Pi 



p1 + ^ 



I4\ In 



Pi 



(3-7) 



Here the second term in the right-hand side of Ea. H3-7() is nothing but the deflection 
angle 5(pi^AB of the mass rm neutrino due to the gravitation. In practical problems. 



it will be sufficient to make approximate 1 + ( — 



2 + 



mjC 



^ ~ and 



PiC J ' \ ^ 

cosh^^ (Rk) — log(2i?fc) in equations (|3-6|) and (|3-7|) . Under these approximations 

2 4 

we can write the defiection angle as 64>i^AB = 2rg(l + ^|^) 

Substituting Ea. ()3-7() for Ea. (|2-9() . the phase function in Ea. p-2|) can be evalu- 



*•* The deflection angle 



:,AB will tend to that of photon 50^'^ ~ 2fg for nii 



0,Rk 
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ated as 



LiCpAB + Sir{B, A) = PilAB + Li5(j)i^AB + Li < -S(j)i^AB log (aRaRb) -fg\. (3-8) 

Here, taking Lifg — ( -f- 2E~) ^9 LiO(j)i^AB — 2rg (— j into account, we finally 

derive the phase difference 

AS^^^{lA,B-rg), (3-9) 

within the approximation up to the order Arn^ = rn^—mf. It may be surprising that 
the AS looked as if it does not depend on Ra, Rb but on the Schwarzschild radius 
rg of the object, around which the neutrinos pass through. In reality, however, these 



parameters are included in the definition of Iab = Ylk=A B y ^\ ~ Po- 
ii) Si for the interior Schwarzschild solution. 



For the neutrino passing through an astrophysical object, we can suppose two 
types of situations. One is the case such that neutrinos pass across the object on 
a long way from ^4 to -B. The other is the long-baseline accelerator experiment on 
the Earth. In both cases, we have to consider the matter effects in addition to grav- 
itational effects. Then, the phase difference AS calculated from gravitational effect 
only should be understood as a supplementary result. 




Fig. 2. The path inside a homogeneous sphere. 



As for the second case, however, it may be worthwhile to show the result obtained 
under the condition that the Earth is a homogeneous sphere of incompressible fluid 
with radius R. The vierbeins are, then, given by 




(3-10) 



The = sin and tg = r are common to the exterior Schwarzschild solution. 
Let us consider the situation of Fig. 2 for the trajectories of neutrinos; that is, Ra = 

*' If we may put <f)AB ~ X^fe^A s then AS can be calculated directly from Ea. H3-6|l . 

The AS in this case is simply obtained by the substitution rg rg{l — log(4_RA^_B)} in Ea. l3-9ll . 
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Rb = R — po ^ rg. Then the fohowing can be obtained (Appendix B): 

AS ~ Iab 1 -] 

In this case, we have to read Eas. H2-6() and ()2-8|) as those of semiclassical approxi- 
mation, since we deal with that the gravity remains on the surface of the Earth. 

§4. Summary and Discussion 

In this paper, we have studied mainly two problems of neutrino oscillation: 
the formulation of neutrino wave function in gravitational fields, its application to 
simple examples of gravitational fields. First, to formulate the wave function in 
gravitational fields, we followed the semiclassical method for the Dirac equation, 
which was discussed by Pauli in a problem of electrons under electromagnetic fields. 
The resultant Dirac wave function is obtained from that of the flat spacetime by 
substituting the phase for Hamilton's principal function of the neutrino in the curved 
spaceitme, providing the observation point is located at the fiat spacetime. We can, 
thus, recognize that the heuristic phases used in curved spacetimes are nothing but 
the result of this approximation. 

Secondly, we attempted to apply our formula to the neutrinos with two gener- 
ations for the gravitational fields in outside of a static uniform sphere and for one 
in inside of the sphere. The latter is applicable, for example, to such an experiment 
as the long-base line, although the effect is examined to be too small in practical 
problems. On the other side, in the former case, the gravitational effect yields the 
probability ()2-8() characterized by the oscillation phase ^ = ^{Iab — fg) with the 
mixing angle 0, where = ^"^^a is the oscillation scale in the vacuum. The result 
says that if the neutrinos started from a source object such as a supernova passes 
through a gravitational area by a star with the scale of sun, then the effective dis- 
tance will be shifted by rg ~ 'ikm. This looks like that the effect of Vg may be ignored 
compared with the distance ~ 10^ light-year, which the neutrinos travel across. We 
must be careful, however, on this problem, since Iab + fg and its remainder divided 
by Iq gives the same oscillation effect; and so, the long distance Iab play the role 
of well mixing of neutrinos between generations. Then, the rg dependence due to 
the pass dependence of neutrinos may give an interesting observational effect in the 
oscillation. 

We note that the gravitational effect in the oscillation phase can be understood 
in terms of the gravitational lensing effect. Indeed, one can write the effect as 
j^rg = ^poAdcj) with AScj) = 64>i — 64>2- Here, 64>i,{i = 1,2) are deflection angles 
for mass eigenstates of neutrinos, which relate to the deflection angle of the photon 
S(f>P^ by dcpi = 6(I)P^{1 + ^■^•i )• If we can observe both deflection angles of a massive 
neutrino and that of photon from the neutrino source object, one can determine Z\m? 
in principle, although such an observation will be very hard in practice. We also note 
that we can introduce the dynamical time tj, the time as a function of r, of each 
particle by ^'^'Ig'^'* = similar to the with deflection angle. The Ea. H3-8() gives 



(3-11) 
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us the difference of these dynamical times as At = tAB2 — tABi = 

in which the gravitational correction is again very small in practical problems. 

Throughout this paper, we discussed the neutrinos under the influence of grav- 
itational fields only provided the mixing angles were constants; that is, that we did 
not take the matter effect into account. We also note that the oscillation phases 
derived in the Schwarzschild spacetime under the condition pQ ^ Vg can not be ap- 
plied for neutrinos passing close to heavy astrophysical objects such as neutron stars, 
though Eas. H3-2|) and 1)3 ■3() hold. In such a case,^ it becomes necessary to develop 
another approximation by taking the spin connections into account. Another basis 
of this paper is the equivalence principle, the universality of gravitational coupling. 
In our approximation, the gravitational effect appears only through Vg. Thus, the 
breakdown of the universality by rg — > rgi simply leads to the substitution A(m?rg) 
for AmPrg, for example, in Ea. ()3-9() . Together with the extension of generations, 
there remain interesting future problems according to this line. 
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Appendix A 

Spin 1/2 particles in curved spacetime 

The curved spacetime can be described by the metric g^i, = rjape^^^ ei^'^ with the 

vierbein {e|r^}, (a, = 0, 1, 2, 3), where r]af3 is the Minkowski metric. The vierbein 
is supposed to satisfy the parallelism condition 

V^e(-) ^ a^eW - r.^^eW + u;(")(p)^e(^) = (A-l) 
; that is, in terms of the inverse vierbein e^^^e^f^ = 6q, 

Vm^h = d,el^ + r^^e^^ - o.(^\„)^e^^) = 0. (A.2) 
Then, the Dirac equation in the curved spacetime has the form 

[in7"ef^)Z)^ - mc]ij = 0, (A-3) 

where D^il) = {d^ — ^a'^^uj(a){f3)^i)4'7 (o""^^ = §[7", 7^]) is the covariant derivative 
acting on the spinor tp. Throughout this paper, we consider the case such that the 
spin connections {w(q)(/3)^} are negligible compared with the other terms in the Dirac 
equation, and so the "0 satisfies 

[i^7"ef„)5^-mc]V' = 0. (A-4) 

To obtain an approximate solution of this equation, let us apply the semi-classical 
method by Pauli,"^ the Dirac counterpart of the WKB approximation; and we put 

iS{x) 

ip = e~i^u, (A-5) 
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where S and u are a scalar phase function and a four components spinor, respectively. 
Substituting this expression for Eq. HA-4|) , we have 



[ih^el-^ {d^ + -d^S) - mc]u = 0. (A-G) 

In non-relativistic Schrodinger equation, the semiclassical approximation is carried 
out by expanding S in powers of j. On the other hand, in the Dirac equation, the 
^ expansion is applied to u instead of that is, 

u = no + ui + U2 + -- - . (A-7) 

Substituting the expansion to (EH, we obtain the successive equations 

(7"ef„)9^5 + mc)no = 0, (A-8) 
(7"ef„)a^S + mc)ui = -^^ef^^a^uo (A-9) 

and so on, at each order of h. Multiplying (7"e|^^^5^5 — mc) from the left-hand side 
of Eq. (|A-8p , we obtain the Hamilton- Jacobi equation 

g^'^d^Sd^S - m^c^ = 0. (A- 10) 

Using the S and a 4-spinor w, the uq can be expressed as 

uoocA^{S)w, (A-11) 
-7"ef ^S^S" + mc 

^-(*' = ■ 

where A_{S) is the projection operator for solutions of Ea. (|A-8| ). 
Now, under the lowest approximation 

li ~ uoe^'^ (A-13) 

Ea.(|A-9j) leads to e['^^5^(uo7°^to) = e1^a)^ti{u7°'u) = by taking uo{j°'e'^^^d^S + 
mc) = into account. On the other hand, since we have disregarded the spin 
connection, the inverse vierbein satisfies V^e^^^ = diJ.G\^a)~^^pi.p^1a) ~ Thus, we can 
write the continuity equation in the covariant form e^^j5^(n7"n) = V^(ne|'^j7"'u) = 
0. This says that the above approximation ()A-13|) is appropriate for the spin 1/2 
particle in the gravitational field under consideration. 

Appendix B 

The phase function inside a homogeneous sphere. 



The phase functions Si have the form H3-2() in this case too; and the counterpart 
of Ea. ()3-3() becomes 



-2 



(B-l) 
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Now, let us consider the case Tg <^ pQ < r < R; then, within the approximation up 
to the Sir can be integrated to the following form: 



Sir — i 



dr 



1 _ '9' 



3f„ 



1 



1-^ + 



2i?, 2_R2 



i?3 



+ 



g2 j,2 



1 - ^ 



2 2 

m,- c 



2L,; 



Vr^ - 1 



COS 



i?2 



+ 



^1 iv^F^ 



+<!+{- 



''^ {ie - 1)1 



(B-2) 



In the practical neutrino experiments connecting two points on the earth, the 
long baseline accelerator experiments, we may put i? ~ 1, since neutrinos pass 
through the ground near surface of the earth. Then the last term in the right-hand 
side of Ea. ()B-2|) can be negligible, and we obtain 



Sir — 2Z/7, 



Vr^-1- cos-1 (4) U f — Y XVr^-1 

RJ } \cpij 2R 



■ PilAB - Li4>AB + 



TglAB I f E 



R Pi \c 



(B-3) 
(B-4) 



where we have identified the geometrical angle 2 cos 



-1 (1 
R 



with the physical angle 



4'ab{= ~^^)' since the spacetime is almost flat on the earth. From this expression 
for Sir, (i = 1, 2), we can derive 



2^3 



AS 



2E 



(B-5) 



Neutrino Oscillations in Gravitational Fields 



11 



References 

1) M. N. Butler, S. Nozawa, R. A. Malaney, and A. I. Boothroyd, Phys. Rev. D 47 (1993) 
2615 

K. lida, H. Minakata and O. Yasuda, Mod. Phys, Lett. A8 (1993) 1037 
A. Halprin and C.N. Leung, Phys. Rev. D53 (1996) 5365 

2) Y. Mellier, Gravitational Lensing, (Springer, 1990) 

P. Schneider, J. Ehlers and E. E. Falco, Gravitational Lenses , (Springer, 1992) 

3) R. Escribano, J. -M. Fere, D. Monderen and V. Van Elewyck, |h'ep-ph/99105'To| v2 28 Oct 
1999 

4) C. Y. Cardall and G. M. Fuller, hep-ph/9610494 Oct 1996 

N. Fornengo, C. Giunti, C.W. Kim and J. Song, Phy. Rev. D56 (1997) 1895 

G. Papini, gr-qc/0110056 (2001) 

Y. Kojima, Mod. Phys. Lett. All, 2965 (1996) 

5) After publishing our paper, we have known the following papers: 

D. V. Ahluwalia and C. Burgard, Gen. Rel. Grav. 28 1161-1170 (1996) 
D. V. Ahluwalia and C. Burgard, Phys. Rev. D57 4724-4727 (1998) 
C. Y. Cardall and G. M. Fuller, Phys. Rev. D 55 (1997) 7960 

R. Escribano, J. -M. Fere, D. Monderen and V. Van Elewyck, Phys. Lett. B512 (2001) 8 
See also: 

K. Konno and M. Kasai, Prog. Theor. Phys. 100:1145-1157(1998) 

6) L. Stodolsky, Gen. Rel. Grav. 11 (1979) 391 
J. Wudka Phys. Rev. D64, 065009 (2001) 

T. Bhattacharya, S. Habbi, and E. Mottola, Phys. Rev. D59 067301 (1999) 

7) W. Pauh, Helv. Phys. Acta 5 (1932), 179: See also 

W. Pauli, General Principle of Quantum Mechanics, §23, Springer- Verlag Berlin Heidel- 
berg New York (1980). 

8) C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation, (W. H. Freeman, 1973) 

9) L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, %101, 4ih Edition 
(Butterworth-Heineman Ltd., 2003) 

10) L. D. Landau and E. M. Lifshitz, Mechanics, §^7, 3rd Edition (Butterworth-Heineman 
Ltd., 1982) 

11) D. Piriz, M. Roy and J. Wudka, |hep-ph/9604403 Apr 1996 



